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ABSTRACT 


The  described  experiment  and  the  evaluation  of  its  results  were 
designed  to  assess  the  significance  of  the  sediment  cracking 
occurring  during  piezometer  and  heater  probe  insertion,  and  the 
influence  of  the  cracking  on  the  subsequent  excess  pore  pressure  and 
thermal  fields  to  be  measured  in  the  In  Situ  Heat  Transfer 
Experiment  (ISHTE).  The  I3HTE  is  part  of  the  Subseabed  Disposal 
Frogram  ( SDP )  managed  by  Sandia  National  Laboratories  Albuquerque 
kSNLA)  and  funded  by  the  Department  of  Energy.  Observed  sediment 
radial  cracking  was  estimated  to  penetrate  to  about  one-half 
penetrator  diameter.  Measured  penetrator  insertion  pore  pressures 
and  excess  pore  pressure  dissipation  indicate  that  sediment  crackincr 
has  no  identifiable  influence  at  the  penetration  depth  of  the 
sensors.  Proper  prediction  of  the  dissipation  of  insertion  pore 
pressures  was  found  to  require  incorporation  of  a  smear  factor, 
accounting  for  sediment  remolding  at  the  piezometer  wail.  An 
analytical  model  to  concurrently  describe  the  rise  in  excess  pore 
pressures  due  to  sediment  heating  and  the  pore  pressure  dissipation 
radially  away  from  the  heat  source  is  described  and  the  results  of 
its  application  to  the  experiment  data  evaluated. 
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1.0  INTRODUCTION 


1 .  1  Purpose 

This  report  presents  the  results  of  a  laboratory  simulation, 
conducted  at  atmospheric  pressure,  of  the  planned  In  Situ  Heat 
Transfer  Experiment  (ISHTE).  The  purposes  of  this  laboratory 
simulation  were  to: 

a)  determine  the  depth  of  sediment  cracking  result ina  from  the 
insertion  of  probes  into  the  sediment, 

b>  measure  and  evaluate  the  excess  pore  pressures  generated 
during  the  insertion  of  probes  into  the  sediment,  and 

c  measure  and  evaluate  the  excess  pore  pressures  generated  by 
the  application  of  a  thermal  gradient  in  the  sediment. 

1*2  Background 

The  Subseabed  Disposal  Program  (SDP)  is  studying  the 
feasibility  of  sequestering  high  level  nuclear  waste  in  fine-grained 
deep-sea  formations  of  the  world  ocean  basins  (Hollister  et  al., 
1^81).  The  ISHTE  field  test  will  provide  data  on  the  thermal,  fluid 
dynamic  and  thermochemical  response  of  in-situ  seabed  sediment  for 
use  m  verification  of  laboratory  experimental  approaches  and 
computer  models  of  waste/'sediment  thermal  interaction.  Durincr  the 
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test,  the  sediment  response  will  be  studied  by  conducting  the 
following  series  of  measurements  (Fercival,  1983): 

a)  Thermal  conductivity; 

b }  Temperature ,  to  develop  a  "picture"  of  the  thermal  field; 

c)  Excess  pore  water  pressure; 

d)  Ion  migration; 

e)  Sampling  of  pore  water  constituents; 

f)  Vane  shear  strength;  and 

g;  Fermeabilities  and  undrained  shear  strengths  of  core 
samples . 

These  tests  will  be  performed  from  the  ISHTE  platform;  both 
and  after  implant  of  an  isotopic  heat  source. 

To  prepare  for  design  of  the  ISHTE  experiment,  the  ISHTE 
Simulation  Experiment  (ISIMU)  was  performed  on  a  0.287:1  scale  model 
of  the  platform.  The  first  test  (ISIMU-I)  was  conducted  from 
November  2  thru  December  15,  1961,  at  the  David  Taylor  Naval  Ship 
Research  and  Development  Center  (Hiller  et  al.,  1982).  In  ISIMU-I, 
the  various  sensors  were  implanted  in  a  re-constituted  iilitic  "red 
clay  sediment  (Silva  et  al.,  1983).  The  sensors  and  sediment 
sample  were  placed  in  a  large  pressure  vessel  and  pressurized  to  55 
MPa,  to  simulate  the  deep  sea  environment,  for  a  one  month  period. 
The  primary  objectives  of  the  experiment  included  (Percival,  1983): 

a;  Determination  of  the  transient  temperature  distribution  in 
the  sediment; 
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b)  Measurements  of  the  thermal  conductivity  of  the  sediment; 

c>  Measurements  of  transient  pore  pressures; 

d  Measurements  of  transient  sediment  response  to  the  applied 
pressure ; 

e;  Measurements  of  the  sediment  shear  strength; 

f )  Studies  of  the  pore  water  chemistrv;  and 

cn  Post-test  studies  of  the  sediment  mineralogy,  chemistry,  and 
structure. 

ISIMU-I  produced  two  unexpected  results.  First,  for  a  given 
heater  power  level,  the  maximum  temperature  of  the  surrounding 
sediment  was  approximately  25  percent  lower  than  that  predicted 
-based  on  pretest  measurements  of  the  sediment  properties  (Percival, 
1983).  The  second  unexpected  result  was  that  the  excess  pore 
pressure  developed  on  insertion  of  the  pore  pressure  probe  in  the 
near-field  (15  mm  from  heating  element)  was  approximately  one  half 
that  of  the  far-field  pore  pressure  probe  (342  mm  from  heating 
element;  (Bennett  et  al.,  in  press).  Cracking  of  the  sediment 
surface  was  noted  to  accompany  the  insertion  of  the  Electric  Heat 
source  ( EHS ) ,  the  pore  pressure  probes,  and  the  geochemistry 
sampler.  These  sediment  cracks  were  hypothesized  to  penetrate  to 
the  vicinity  of  the  heater  source  and  the  pore  pressure  sensors  and 
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to  be  responsible  for  the  unexpected  temperature  and  excess  pore 
pressure  dissipation  results.  ISIMU-I  also  produced  information  on 
the  impact  of  sediment  heating  on  the  excess  pore  pressures.  Upon 
heating,  the  near-field  piezometer  displayed  a  very  short  term  lag 
m  response  t<4  min)  followed  by  a  short  term  C  <40  min)  increase  in 
pressure  with  a  subsequent  long-term  (*90  hrs)  exponential  decrease 
m  pressure.  The  far-field  piezometer  showed  a  significant  short 
term  lag  08  hrs)  followed  by  a  long  term  0120  hr)  pressure 
increase.  The  two  pressures  for  the  very  long  term  (>  120  hrs) 
leveled  off  at  approximately  the  same  values. 

ISIMU-I  then,  while  providing  answers  to  many  questions,  also 
left  questions  unanswered  and  raised  new  ones.  This  report 
describes  a  follow-up  experiment  designed  to  help  resolve  some  of 
the  questions  remaining  after  ISIMU-I. 

1.3  Approach 

To  further  evaluate  the  ISIMU-I  results,  an  additional 
simulation  experiment  (ISIMU-II)  was  performed.  This  test  was 
conducted  from  April  2  to  April  5,  1984,  at  Sandia  National 
Laboratories  in  Albuquerque  (SNLA).  Since  only  the  thermal  and  pore 
pressure  response  characteristics  of  the  sediment  were  of  interest, 
ISIMU-II  did  not  involve  pressurization  to  the  deep-sea  environment 
or  the  same  level  of  instrumentation  as  ISIMU-I.  Two  pore  pressure 
probes  were  placed  within  15  mm  of  a  heating  probe  in  a  0.46  m 
diameter  tank  filled  with  re-constituted  illite  sediment  (the  same 
sediment  used  in  ISIMU-I).  The  heater  probe  and  the  tank  walls  were 
instrumented  with  thermocouples.  Pore  pressures  were  monitored 
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during'  the  individual  probe  insertions ,  and  during  subsequent 
heating . 

This  report,  describing  ISIMU-II  and  its  results,  is  arranged 
in  four  sections.  The  first  section  describes  the  test  arrangement 
and  data  reduction.  The  next  section  discusses  the  development  of 
models  to  predict:  (1)  the  excess  pore  water  pressure  dissipation 
following  probe  insertion  and  (2)  the  excess  pore  pressure  build-up 
due  to  thermal  effects.  The  third  section  presents  the  results  of 
the  experiment  and  the  application  of  the  proposed  prediction 
models.  The  last  section  summarizes  the  study. 

2.0  TEST  ARRANGEMENT 

2.1  Sediment  Sample  Preparation 

The  sediment  sample  used  in  the  ISIMU-II  experiment  was 
prepared  at  SNLA .  The  sediment  used  was  the  same  pelagic  clay  from 
MF'G-I,  of  primarily  illitic  composition,  used  in  the  I3IMU-I 
experiment.  This  material  had  been  shipped  from  NSRDC ,  Annapolis. 
MD,  to  SNLA  in  55-gallon  drums  and  stored.  This  material  was 
reconstituted  in  a  plaster  mixer  into  a  thick  slurrv  state  and 
poured  into  the  sediment  tank  in  0.15  m  layers.  Vacuum  was  applied 
after  each  layer  had  been  placed  to  de-air  the  material.  The 
sediment  tank  used  in  ISIMU— II  was  a  glass  cylinder  with  a  sand 
filter  at  the  bottom  (Figure  1),  burlap  fabric  filters  on  the  sides, 
and  a  fiber  filter  sheeting  between  the  sediment  and  the  burlap  to 
limit  clogging  of  the  burlap.  The  slurry  sample  was  loaded  in 
increments,  using  deadweights,  to  duplicate  the  consolidation  state 
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achieved  in  ISIMU  I.  The  resulting  sample  was  0.43  m  (17  in)  in 
diameter  and  0.71  m  (28  in)  in  height. 

2.2  Condition  of  Sample  at  Testincr 

Cores  were  taken  at  the  conclusion  of  the  test  and  water 
contents  and  miniature  vane  shear  strengths  measured  on  those  cores 
The  results  of  these  tests  on  the  cores  are  reported  and  discussed 
later  in  this  report.  Tests  on  the  peripheral  cores  show  the 
sediment  sample  to  be  of  fairly  uniform  water  content  (average  88.0 
percent;  and  undrained  shear  strength  (average  4.1  kPa)  with  depth. 

samples  and  tests  were  not  conducted  to  provide  for 
evaluating  the  variability  in  water  content  and  strength  across  the 
diameter  of  the  test  tank;  however ,  prior  experience  with  this 
sample  preparation  technique  for  ISIMU-I  suggests  this  variation 
will  not  be  significant  (Silva  et  al . ,  1983). 

2.3  Pore  Pressure  Probes 

The  pore  pressure  probes  are  of  the  same  design  as  those  used 
in  ISIMU-I  and  to  be  used  on  the  ISHTE  platform.  The  probe  sensing 
element  is  mounted  on  a  7.94  mm  (5/16  in.)  diameter  titanium  tube 
about  1  m  in  length.  The  probe  tip  is  a  cone  100  mm  in  length  and 
7.9  mm  diameter  at  its  base  (5.3  degree  cone  angle),  designed  to 
induce  two-dimensional,  or  lateral,  deformation  during  penetration 
into  the  sediments  (Bennett  et  al.,  1981  and  in  press).  The  excess 
pore  pressures  are  sensed  at  a  13  mm  long  (1/2  in.)  cylindrical 
porous  corundum  stone  located  160  mm  from  the  probe  tip.  The  pore 
pressures  sensed  at  this  porous  stone  are  transmitted  through  the 
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titanium  tube  to  a  variable  reluctance  differential  pressure 
transducer  located  1  m  above  the  porous  stone  (Bennett  et  al.,  1980, 
1983,  and  in  press). 

2.4  Insertion  Experiment 

ISIMU-II  was  intended  to  duplicate  the  heater  probe  and  near¬ 
field  pore  pressure  probe  insertions  performed  in  ISIMU-I.  In 
ISIMU-I ,  the  near-field  pore  pressure  probe  was  positioned  15  mm, 
wall-to-wall,  from  the  heater  probe,  and  it  had  been  planned  to 
position  the  porous  stone  at  the  same  elevation  as  the  center  of  the 
heating  element.  However,  post-test  inspection  revealed  that  the 
nea.r-field  piezometer  had  not  attained  the  planned  penetration,  but 
had  stopped  101  mm  short  {Miller  et  al . ,  1982).  Therefore,  in  order 
to  duplicate  the  original  test  conditions,  the  porous  stones  of  the 
probes  used  in  ISIMU-II  were  both  positioned  101  mm  above  the  mid¬ 
plane  of  the  heat  source  (Figure  2). 

The  original  plan  of  ISIMU-II  called  for  installation  of  the 
heater  probe  at  the  center  of  the  tank.  A  100  mm  (4-in.)  thick 
block  of  Teflon,  with  3  bored  guide  holes  for  the  heater  and  two 
pore  pressure  probes,  was  fixed  over  the  top  of  the  tank.  The  first 
piezometer  was  saturated  by  first  immersing  the  probe  tip  including 
porous  stone  in  the  150  mm  (5-3/4  in.)  of  water  overlying  the 
sediment  surface.  A  partial  vacuum  was  then  applied  to  the  positive 
bleed  port  of  the  transducer,  standing  some  0.7  m  above  the  water 
surface,  and  water  thus  drawn  up  into  the  transducer  to  saturate  the 
system.  On  installation,  this  first  probe  exhibited  a  response  time 
of  3  minutes,  which  was  considered  to  be  too  long  for  a  properly 
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saturated  pore  pressure  probe.  This  probe  was  removed  the  following 
day  and  replaced  with  a  dummy,  pointed,  7.94  mm  (5/16-in.)  diameter 
stainless  steel  rod,  to  fill  the  hole. 

To  minimize  the  disturbance  and  the  discontinuity  of  the  dummy 
rod,  the  Teflon  guide  block  was  reoriented  90  degrees  and  the 
proposed  location  of  the  heater  moved  44  mm  (1-3/4  in.)  off  tank 
center.  For  the  first  pore  pressure  probe,  denoted  as  P-3,  the 
saturated  porous  stone  and  cone  tip  were  attached  under  the  water 
suriace  to  facilitate  system  de-airing.  A  saturation  sleeve  was 
installed  over  the  porous  stone  and  water  was  driven  through  the 
sleeve  and  stone  and  up  the  probe  tube/shaft  to  fill / saturate  the 
transducer  positive  side.  This  probe  was  then  inserted  over  a 
period  of  8  seconds  and  exhibited  a  response  time  of  40  seconds  to 
peak  insertion  pressure,  which  was  considered  indicative  of  a 
satisfactorily  saturated  system.  The  excess  pore  pressures 
generated  by  this  pore  pressure  probe  were  allowed  to  dissipate  for 
156  minutes  (analysis  of  the  data  shows  that  95%  dissipation 
occurred  at  about  90  minutes).  The  heater  (23.7  mm  diameter)  was 
then  inserted,  again  guided  by  the  Teflon  block,  with  the  heater 
wall  15  mm  from  the  pore  pressure  probe  wall  (Figure  2).  Insertion 
pressures  generated  by  the  unpowered  heater  probe  were  allowed  to 
dissipate  over  19  hours.  The  second  piezometer  probe,  designated  p- 
2,  was  saturated  in  the  same  manner  as  P-3,  was  installed  on  the 
opposite  side  of  the  heater  from  P-3  (Figure  2),  also  exhibited  a  40 
sec  response  time  to  peak  insertion  pressure,  and  was  allowed  118 
minutes  for  excess  pore  pressure  dissipation. 
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2.5 


Heater  Experiment 


Seven  thermocouples  were  located  in  the  sediment  filled  tank  in 
the  positions  noted  in  Figure  2.  The  number  and  location  of  each 
are  as  follows: 


0. 

1. 

2. 

3. 

4. 

5. 

b  . 


At  the  tank  wall  in  the  sediment,  280  mm  below  mudline. 
On  the  heater  element,  50  mm  below  mid-plane. 

On  the  heater  element,  at  mid-plane. 

On  the  heater  element,  50  ram  above  raid-plane. 

On  the  heater  element,  at  mid-plane. 

At  the  tank  surface,  in  air,  at  ambient  temperature. 

At  the  tank  wall  in  the  sediment,  280  mm  below  mudline. 


After  the  excess  pore  pressures  from  the  insertion  of 
piezometer  probe  P-2  had  dissipated,  the  heater  probe  element  was 
energized  to  about  29.4  watts.  This  power  input  was  maintained  for 
251  minutes  while  developing/dissipating  excess  pore  pressures  were 
monitored  on  pore  pressure  probes  P-3  and  P-2.  The  heater  response 
portion  of  the  experiment  was  terminated  after  251  minutes  due  to  a 
power  failure.  Even  though  the  power  was  restored  within  5  minutes, 
the  hiatus  in  the  energized  condition  was  sufficient  to  adversely 
impact  the  analysis  of  subsequent  data  (given  the  existing 
analytical  models)  therefore  continuing  the  experiment  after  the 
power  failure  was  without  merit.  However,  sufficient  data  were 
collected  during  the  251  minutes,  before  the  power  failure,  to  make 
the  heater  response  portion  of  the  experiment  meaningful. 
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2.6  Data  Reduction 


The  test  data  were  recorded  using  a  Fluke  data  logger  and  a 
Hewlett-Packard  strip  chart  recorder.  The  raw  data  consisted  of  the 
time  (days,  hours,  minutes,  and  seconds)  at  which  the  measurements 
were  made,  voltage  output  for  each  transducer  and  temperature 
indicated  by  each  thermocouple.  A  tabulation  of  the  data  is  given 
in  Addendum  I. 

The  data  were  stored  on  a  floppy  disc  and  reduced  using  a 
Hewlett-Packard  model  85  computer.  The  times  were  adjusted  to 
reflect  the  total  elapsed  time  since  the  beginning  of  the  test.  The 
voltages  were  transformed  to  pressures  using  the  following 
equations: 

for  probe  P-3, 

Pressure  (psi)  =  0.0111  +  3.9982  a  (volts  +  0.1505)  ...  (l) 

and  for  probe  P-2, 

Pressure  (psi)  =  0.0391  +  3.9950  A  (volts  +  0.1071)  .  .  .  (2) 

The  values  0.1505  and  0.10/1  volts  represent  the  zero  offset 
correction  for  the  probes  inserted  in  the  sediment. 

The  variation  in  pressure  for  the  two  probes  over  the  duration 
of  the  test  are  shown  in  Figure  3.  The  pressures  for  the  individual 
test  events  are  shown  in  Figures  4  to  7.  It  should  be  noted  that 
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the  small  peak  at  the  end  of  Figure  4  is  an  actual  event.  A  lucite 
collar  on  the  pore  pressure  probe  (P-3)  had  to  be  removed  in  order 
to  provide  clearance  for  insertion  of  the  heater  probe.  Durincr 
removal  of  this  collar  the  probe  P-3  was  twisted  gently  and  only 
slightly.  Disturbance  of  the  sediment  during  this  slight  movement 
resulted  in  a  0.31  kPa  (0.045  psi)  excess  pore  pressure.  The 
beginning  of  Event  2,  the  heater  probe  insertion,  was  postponed 
about  18  minutes  to  allow  for  dissipation  of  the  excess  pore 
pressure  generated  by  the  disturbance. 

In  addition,  it  should  be  noted  that  the  initial  pore  pressure 
drop  shown  in  Figure  7  is  a  result  of  the  time  lag  between  the 
initiation  of  heating  and  the  pore  pressure  response  at  the  probe  15 
mm  away.  The  pressure  drop  is  a  result  of  pore  pressure  dissipation 
carried  over  from  Event  3  (second  pore  pressure  probe  insertion).  A 
complete  listing  of  the  reduced  data  is  given  in  Addendum  II. 

3.0  GENERAL  THEORY  AND  MODEL  DEVELOPMENT 

3.1  Pore  Pressure  Dissipation 

Several  theories  (Acar  et  al . ,  1982;  Baligh  and  Levadoux,  1980; 
Levadoux  and  Baligh,  1980;  Randolph  et  al.,  1979;  Soderberg,  1962; 
Torstensson,  1977;  Tumay  and  Acar,  1984,  and  Wroth  et  al . ,  1979) 
have  been  advanced  pertaining  to  the  dissipation  of  excess  pore 
water  pressures  induced  by  the  expansion  of  a  cylindrical  cavity. 

The  governing  differential  equation  in  two  dimensions  is  as  follows: 


dU  =  c. 

l 

a 

i - 

31 

*T)| 

at  h 

r 

ar 

\  dr  / 
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where , 


u  =  excess  pore  pressure, 
t  =  time, 

r  =  radial  distance,  and 

ch  -  horizontal  coefficient  of  consolidation. 

Due  to  the  usual  lack  of  horizontal  coefficient  of 
consolidation  data,  ch,  the  constant  ch  is  usually  replaced  by  the 
coefficient  of  vertical  consolidation,  c v;  which  is  normally  equal 

to  or  less  than  c^  (Soderberg,  1962).  The  coefficient  of  vertical 
consolidation  is  defined  as: 

k  ( 1  +  e ) 

v  ~a - « .  . . 

v  1 

where , 

k  =  permeability 
e  =  void  ratio, 

av  =  coefficient  of  vertical  compressibility,  and 
y  =  unit  weight  of  water 

The  most  commonly  used  method  of  solving  Equation  (3)  involves 
expressing  the  excess  pore  pressure  in  terms  of  a  dimensionless 
dissipation  factor  U  (u/u^^)  which  then  becomes  a  function  of  a 
dimensionless  time  factor  expressed  as. 
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where , 


radius  of  the  embedded  cylinder. 


The  functional  relationship  between  U  and  T  depends  primarily  on 
how  one  characterizes  the  sediment.  Past  pore  pressure  probe 
studies  (e.g.,  Bennett  et  al.r  1984)  have  relied  on  the  solutions  of 
Soderberg  (1962)  and  Randolph  et  al.,  (1979),  which  were  developed 
for  a  foundation  pile  embedded  in  either:  (a)  an  elasto-plastic 
medium  or  (b)  a  viscous  medium.  Using  the  general  form  of  the 
solution  one  can  obtain  Tr  values  for  various  levels  of  dissipation. 
The  most  common  procedure  uses  the  fifty  percent  dissipation  level 
which  yields  a  T  approximately  equal  to  one.  Then  Equation  (5)  can 
be  rewritten  as. 


Ch  = 


r  5  o 


v  6 ) 


5  0 


where  t  represents  the  actual  time  at  which  fifty  percent 
dissipation  has  occurred,  as  obtained  from  a  U  versus  locr  t  plot. 

With  c^  thus  determined  from  Equation  (6),  U  can  be  defined  in 
terms  of  an  exponentially  decreasing  function  of  T  .  This  method 
normally  employs  a  series  (Scott,  1963);  however,  for  the  sake  of 
simplicity  only  the  first  term  in  the  series  will  be  used.  In 
addition,  the  expression  is  written  in  terms  of  real  time  as, 

~3t 


u 


u 


=  e 


v  / 1 


max 
where , 
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The  value  of  m  depends  upon  the  assumptions  used  in  the  theory. 
Scott  (1963)  introduced  the  concept  of  a  "smear"  factor  for  which 
the  parameter  m  is  used  to  reflect  the  reduction  in  permeability  at 
the  boundary  due  to  remolding.  This  concept  can  be  used  to  account 
for  the  disturbance  between  the  probe  and  the  surrounding  sediment. 

3.2  Thermally  Induced  Pore  Pressures 


An  induced  thermal  field  can  have  a  pronounced  effect  on  the 
pore  pressure  response  of  a  sediment.  If  the  temperature  increase 
is  rapid,  a  significant  excess  positive  pore  pressure  can  develop, 
due  to  the  greater  volumetric  expansion  of  the  sediment  pore  water 
than  that  of  the  mineral  solids.  The  thermally  induced  pore 
pressure  change  for  the  undrained  condition  can  be  expressed  in 
terms  of  the  following  equation  developed  by  Mitchell  (1976): 


Au 


where , 


AT 


m 
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n  ( * 
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«  ) 
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DC 

St 


(9) 


AT  =  temperature  change, 

~  coefficient  of  compressibility  of  the  soil  structure 
from  the  rebound-reloading  curve  of  a  one-dimensional 
consolidation  test, 

n  =  porosity, 

ocs  =  thermal  coefficient  of  volumetric  expansion  for  the 
mineral  solids, 

ocw  -  thermal  coefficient  of  expansion  of  the  pore  water,  and 

ocst  =  thermal  coefficient  of  expansion  for  the  soil  structure. 


In  considering  the  problem  of  the  pore  pressures  at  the  pore 
pressure  probe,  one  can  think  of  the  transient  thermal  flux,  i.e., 
when  the  temperature  "front"  first  contacts  the  probe,  as  being  of 
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the  undrained  case.  However,  in  the  long  term  case,  as  temperatures 
attain  the  "steady  state"  condition  significant  spatial  temperature 
gradients  form  in  the  sediment  medium  (Percival,  1983).  These 
gradients  result  in  a  dissipation  of  the  excess  pore  water  pressure 
in  the  direction  of  decreasing  temperature,  i.e.,  to  zones  of  lower 
pore  water  pressures.  Therefore,  in  developing  a  model  for  the 
thermally  induced  pore  pressures  at  the  probe  one  must  consider  not 
only  the  gain  in  pore  pressure  with  respect  to  temperature  and  time, 
but  also  the  loss.  In  lieu  of  a  more  rigorous  dissipation  model 
involving  temperature,  the  proposed  model  will  involve  a 
superposition  of  the  pore  pressure  gain  and  loss,  viz.; 


u 


u 


gain 


u 


loss 


1 10) 


where  u  and  u. 

gam  1 

respectively.  The 
function  of  time, 
time  sequence: 


oss  are  represented  by  Equations  (9;  ana  (7), 

physical  sicrnif icance  of'  u  and  u,  ,  as  a 

gain  loss 

can  be  observed  from  Figure  8  and  the  following 


t  =  1 :  u  =  Au 

gain  i 

u,  =  Au,  ~  hu  e  ~  ti) 

loss  i  i 


t  =  2 :  u  .  =  Au  +  Au 

gam  1  2 


u 


=  (Au  -  Au  e  ^(t2  ti))  +  (Au  -  Au  e" 


loss  "  '“ui  *  1  7  T  v“u2  ~  ^u2  e  P(tz  t2)) 


t  =  3 :  u  .  =  Au  +  Au  +  Au 

gam  123 


u 


loss 


=  (Au  -  Au  e  ^(t3  ti);+  (Au  -  Au  e  ^(t3  t2) 

11  2  2 

+  (Au  -  Au  e  ^lt3  3  * ) 

3  3 
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hence  the  following  series  for  any  time,  k: 


u 


gam 


k 

y 

1  -  ji-11 

i  =  l 
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1 11a) 


k 


u 


loss 


( lib) 


i  =  l 


i=l 


substituting  Equations  (11a  and  lib)  into  Equation  (10)  and 
cancelling  the  like  terms  yields  the  total  excess  pore  pressure  at 

time  t,  : 
k 


k 


( 12) 


i=l 


It  should  be  noted  that  in  the  above  formulation,  Au  was  used 
in  place  of  '-ijnax  in  Equation  (7).  This  substitution  is  physically 
correct  in  that  Au  represents  the  maximum  pore  pressure  beincr 
dissipated  for  a  specific  time  increment.  Thus,  Equation  (12)  can 
account  for  thermally  induced  pore  pressures,  as  well  as  the  time 
dependent  dissipation  of  these  pore  pressures. 

Note  here  that  the  coefficient  (3  is  not  the  same  as  the 
coefficient  8  defined  in  Equation  (8).  The  coefficient  8  describes 
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the  excess  pore  pressure  dissipation  at  the  surface  of  the 
piezometer  probe  and/or  heater  probe:  j3  describes  dissipation 
radially  away  from  the  cylindrical  surface  of  radius  r  .  However, 
during  Event  4,  that  period  with  the  heater  turned  on,  the  pore 
pressure  field  radially  around  the  piezometer  probes  is  quite  non- 
symmetric  and  probably  not  described  by  Equation  (8).  The  form  of 
the  dissipation  function.  Equation  (7),  was  assumed  to  remain 
reasonable,  and  a  new  coefficient,  (3,  adopted  to  describe  the 
dissipation  performance  at  the  piezometer  resulting  from  heater 
heating.  A  function  for  (3  has  not  been  identified. 

4.0  TEST  RESULTS  AND  MODEL  VERIFICATION 

4.1  Sediment  Cracking 

One  purpose  of  the  ISIMU-II  experiment  was  to  estimate  the 
depth  of  sediment  cracking  resulting  from  probe  insertion  and  to 
assess  the  impact  of  this  cracking  on  the  excess  pore  pressure 
dissipation  and  on  the  thermal  field.  Insertion  of  the  pore 
pressure  probes  and  the  heater  probe  did  result  in  cracking  of  the 
sediment  surface.  Radial  cracks  occurred  at  30-60°  spacing  around 
each  probe  extending  outward  to  about  2-1/2  probe  diameters.  These 
cracks  appeared  to  extend  to  a  sediment  depth  of  about  one-half 
probe  diameter.  The  crack  pattern  suggests  an  outward  displacement 
of  shallow,  pie-shaped  wedges  of  sediment  during  the  initial  stages 
of  penetration  followed  by  a  less-obvious  heave  of  the  local  area 
during  continued  penetration.  The  shallowness  of  the  crack  system 
suggests  that  it  will  have  a  negligible  influence  on  the  excess  pore 
pressure  dissipation  and  temperature  gradient  in  the  vicinity  of  the 
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heat  source  and  the  pore  pressure  sensors.  The  similar  dissipation 
responses  of  the  two  piezometers,  one  installed  before  the  heater, 
the  other  after,  demonstrates  that  cracking  due  to  heater 
installation  has  not  measureably  influenced  the  response  of  the 
second  piezometer,  P-2. 

4.2  Pore  Pressure  Dissipation 

This  section  will  deal  with  application  of  the  theory  developed 
in  Section  3.1  to  data  presented  in  Section  2.6.  In  order  to 

determine  the  horizontal  coefficient  of  consolidation,  c  ,  iM  = 

h 

necessary  to  plot  the  data  shown  in  Figures  4,  5,  and  6  on  a  semi- 
log  plot  to  evaluate  t10Q  and  then  t5Q  tc.f.  Equation  (6)).  The 
excess  pore  pressure  data  were  first  normalized  by  dividing  bv  the 
maximum  observed  pore  pressure  value.  These  data  are  plotted 

log  time  as  shown  in  Figures  9  and  10  for  pore  pressure 
probes  P-3  and  P-2,  respectively.  The  tests  were  interrupted  at  95 
to  98  percent  dissipation;  therefore,  a  lower,  asymtotic  value, 
i.e.,  complete  dissipation  of  excess  pore  water  pressures,  was  not 
attained.  Thus,  a  lower  bound  value  was  interpreted  and  is  shown  as 
the  coarse  dashed  lines  in  Ficrures  9  and  10.  The  value  of  t 

50 

corresponds  to  the  time  at  which  u/u  is  slightly  crreater  than 
0.5.  An  exact  value  of  0.5  was  not  used  since  the  lower  asymptote 
did  not  reach  0.0.  The  t,_Q  times  are  8.3  and  11.4  minutes  for 
probes  P-3  and  P-2,  respectively.  The  ch  values  obtained  from 
Equation  (6)  using  a  probe  radius  of  0.4  cm  are  3.21  x  10-4  and  2.34 
, . -4  2  . 

x  10  cm  /sec  for  the  respective  probes. 
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Figure  11  is  the  same  type  of  plot,  with  the  exception  that  it 
is  for  the  excess  pore  pressure  dissipation  at  probe  P-3  due  to  the 
insertion  of  the  heater  probe.  The  t^^  time  for  this  event, 
measured  at  probe  P-3,  is  53  minutes  as  compared  to  t^  times  of  8 
and  11  minutes  for  the  probe  insertions.  The  longer  dissipation 
time  for  the  heater  probe  insertion  is  due  primarily  to  the  larger 
cross-sectional  area  of  the  heater  probe  as  compared  to  the  area  of 
the  pore  pressure  probes,  with  dissipation  time  at  the  probe  surface 
being  directly  proportional  to  the  probe  radius  squared  (from 
Equation  (6)).  Further,  the  t^  measured  at  the  pore  pressure  probe 
is  somewhat  longer  than  the  t ^  that  would  be  measured  at  the  heater 
probe  surface  because  the  dissipation  rate  decreases  as  one  moves 
away  from  the  probe  surface.  An  additional  point  of  interest  is  the 
prediction  that  the  pore  pressure  decreases  inversely  with  the 
radial  distance  away  from  the  probe  (Soderberg,  1962),  viz.; 


where 

=  excess  pore  pressure  in  the  sediment  at  a  radial  distance 
r.  from  the  probe, 

u  =  excess  pore  pressure  at  the  probe  of  radius  r 
u  o 

Assuming  the  maximum  excess  pore  pressure  generated  during  heater- 
probe  insertion  to  be  the  same  as  that  resulting  from  the  insertion 
of  probe  P-3  (11.6  kPa  or  1.69  psi),  one  obtains  a  predicted  value 
of  5.2  kPa  (0.75  psi)  at  probe  P-3  due  to  heater  insertion  (using 
rQ  =  12  and  ri  =  15  mm).  (The  above  statement  assumes  that  the 
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change  in  tip  shape,  very  sharp  (pore  pressure  probe)  versus  blunt 
(heater  probe)  and  the  change  in  probe  diameter  from  7.9  mm  to  23.7 
mm,  have  no  influence  on  the  excess  pore  pressure  generated  during 
insertion.  This  predicted  value  compares  favorably  with  the 
observed  value  of  5.8  kPa  (0.84  psi).  The  predicted  pressure 
response  at  probe  P-3  due  to  the  insertion  of  P-2  is  0.69  kPa  (0.10 
psi)  (using  rQ  =  4  and  r..  =  54  mm)  as  compared  to  the  observed  value 
of  0.97  kPa  (0.14  psi).  This  comparison  ignores  the  influence  of 
the  heater  probe,  located  between  the  pore  pressure  probes,  on  the 
generated  excess  pore  pressure  at  P-3. 

The  above  values  are  used  in  the  dissipation  model  (Equation 
(7),  assuming  m  =  1)  to  predict  the  dissipation  rates  shown  in 
Figures  12  and  13  for  probes  P-3  and  P-2,  respectively.  Also  shown 
in  these  figures  is  the  least  squares  exponential  curve  fit  to  the 
data  ( "*  represents  data  point).  The  fit  is  less  than  ideal, 
because  the  intercept  value  is  forced  to  be  unity,  leavinq  only  one 
parameter  i the  slope)  to  be  evaluated.  The  important  point  to  note 
in  Figures  12  and  13  is  the  very  close  agreement  between  the 
predicted  and  observed  values  for  the  short  term  (<  5  minutes) 
behavior.  However,  in  terms  of  medium  to  long-term  behavior  the  two 
curves  differ  considerably,  with  the  predicted  curve  showing  a  much 
dissipation  rate.  Between  the  predicted  and  least  squares 
curves  the  latter  will  be  considered  the  more  correct  in  predictinq 
the  overall  dissipation  behavior.  Therefore,  it  becomes  necessary 
to  "slow  down"  the  dissipation  rate  as  calculated  from  Equation  (7). 
This  can  be  achieved  through  the  use  of  the  m  value  in  Equation  (8). 
Figures  14  and  15  show  the  effects  of  the  m  parameter  in  shifting 
the  predicted  curves  for  probes  P-3  and  P-2,  respectively.  The 
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corresponding  m  values  which  resulted  in  curves  closest  to  the  least 
squares  curves  are  4.6  and  2.3  for  the  respective  probes.  Thus,  of 
the  two  probes,  the  dissipation  curve  for  probe  P-3  is  affected  to  a 
greater  extent  (higher  m  value)  due  to  insertion.  This  finding  is 
consistent  since  P-3  was  inserted  in  "virgin"  sediment,  whereas,  P-2 
was  inserted  in  sediment  disturbed  by  the  prior  insertion  of  the 
heater  probe. 

4.3  Thermally  Induced  Pore  Pressures 

The  thermally  induced  excess  pore  pressures  at  probes  P-2  and 
P-2  were  previously  shown  in  Figure  7.  The  salient  features  to  be 
noted  from  this  figure  are  as  follows: 

1)  the  initial  4  minute  lag  in  the  pore  pressure  response, 

2)  the  significant  pore  pressure  decay  for  times  greater  than 

75  to  100  minutes, 

3)  the  reduced  pore  pressure  response  of  probe  P-2,  and 

4)  the  hyperbolic  shape  of  the  initial  portion  of  the  curves. 

The  first  factor  can  be  directly  attributed  to  the  thermal 
diffusivity  of  the  sediment.  The  same  4  minute  response  lacr  was 
observed  during  the  ISIMU-I  simulation  (see  Section  1.2).  Shifting 
the  curves  toward  the  origin  compensates  for  this  factor.  The 
second  and  third  factors  are  a  result  of  the  previous  insertion 
event.  Figure  3  clearly  shows  that  pore  pressure  dissipation  during 
this  event  had  not  attained  the  asymptotic  minimum  prior  to  heater 
"turn-on".  Thus,  the  previous  dissipation  history  must  be  taken 
into  account.  This  is  achieved  by  fitting  a  least  squares 
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exponential  curve  to  the  "tail"  portion  of  the  curves  shown  m 
Figure  6.  The  correction  is  added  to  the  original  data  as  is  shown 
in  Figures  16  and  17  for  probes  P-3  and  P-2,  respectively.  The 
fourth  factor  is  a  result  of  the  rate  of  temperature  increase  as 
described  by  Equation  (9).  Equation  (9)  predicts  the  rate  of  pore 
pressure  build-up  to  be  directly  related  to  the  rate  of  temperature 
change.  The  temperature- time  curves,  as  recorded  by  thermocouples 
located  on  the  heater  probe,  are  shown  in  Figure  18.  One  observes 
from  this  figure  the  hyperbolic  shape  of  the  curves.  The  actual 
temperature-time  relationship  at  the  piezometer  probes  was  not 
recorded.  However,  the  thermocouples  at  the  tank  wall,  200  mm  from 
the  probes,  recorded  no  significant  temperature  change.  Therefore, 
the  actual  heating  rate  must  lie  somewhere  between  that  shown  in 
Figure  18  and  zero.  Since  the  actual  curve  is  not  known,  the 
hyperbolic  form  will  be  used  with  assumed  parameters. 

The  hyperbolic  relationship  is  used  to  calculate  the  chancre  in 
temperature  for  a  unit  time  step  which  is  then  substituted  into 
Equation  (9)  to  determine  the  rate  of  pore  pressure  increase.  As 
mentioned  in  Section  3.2,  the  pore  pressure  increase  is  affected  by 
dissipation  into  the  surrounding  sediment,  therefore,  the  unit  pore 
pressure  increase  must  be  substituted  into  Equation  (12)  and  summed 
over  the  entire  record  to  obtain  the  actual  pore  pressure  response. 
This  response  is  shown  in  Figures  19  and  20  for  probes  P-3  and  P-2, 
respectively.  In  these  figures  an  asterisk,  "A",  represents 
corrected  data,  as  described  previously,  "  +  "  represents  the 
thermally  induced  pore  pressures  as  obtained  from  the  assumed 
temperature-time  relationship  and  Equation  (9),  and  "0"  represents 
the  predicted  pore  pressure  build-up  using  the  thermally  induced 
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pore  pressure  data  and  Equation  (12).  For  this  calculation,  as  a 
first  approximation,  the  coefficient  {3  was  assumed  equal  in  value  to 
the  coefficient  3  from  Equation  (8).  The  fit  of  prediction  to 
measured  data  shown  in  Figures  19  and  20  was  achieved  without 
further  adjustment  to  the  value  of  f3.  The  computer  program  "MENUDO" 
used  to  create  these  figures  is  given  in  Addendum  III.  The  program 
listing  gives  the  typical  values  for  the  sediment  structural 
properties,  as  well  as,  the  thermal  properties  of  illite  (Bennett  et 
al.,  1984;  Mitchell,  1976;  and  Silva  et  al.,  1982). 

The  authors  note  that  much  work  remains  before  the  temperature- 
pore  pressure  response  model  described  herein  is  accepted  as  a 
reliable,  effective  model.  The  present  model  uses  a  number  of  soil 
behavior  parameters,  e.g.,  c^,  c^,  m^,  at  stress  and  temperature 
levels  where  these  behavior  parameters  are  only  vaguely  understood. 
Back-calculations  of  the  sediment  temperature  at  the  piezometer 
yield  calculated  temperatures  ranging  from  one  order  of  magnitude 
higher  to  one  order  of  magnitude  lower  than  the  measured 
temperatures  at  the  heater  surface,  all  for  reasonable  assumptions 
for  input  parameter  values.  Continued  effort  is  required  to  improve 
our  understanding  of  the  temperature-excess  pore  pressure 
generation-dissipation  phenomena  to  the  working  level. 

4.4  Sample  Evaluation 

After  completion  of  the  test  program.  Dr.  Les  Shepherd,  of 
Sandia  National  Laboratories,  obtained  push  cores  of  the  tank 
sediment  near  the  periphery  of  the  tank  in  sediment  believed  least 
disturbed  during  the  experiment .  These  peripheral  core  tubes  were 
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inserted  before  any  of  the  test  probes  were  withdrawn  to  minimize 
possible  lateral  deformations  and  disturbance  of  the  sediment  during 
coring.  The  peripheral  cores  and  probes  were  then  withdrawn  and  the 
area  of  the  heater  probe  was  then  over-cored.  Two  of  the  peripheral 
cores  and  the  heater  overcore  were  then  tested  at  SNL,  and  one 

peripheral  core  transported  as  handcarried  luggage,  was  tested  at 
NORDA. 

Water  contents  measured  on  two  cores  tested  at  SNL  (without 
transportation)  show  the  water  content  to  be  reasonably  constant 
with  depth,  varying  from  85  to  90  percent  with  an  average  for  the 
two  cores  of  87.6  percent  (Figure  21).  Water  contents  on  the  one 
core  tested  at  NORDA,  after  transport  by  plane  and  auto,  were 
essentially  the  same  ranging  from  87  to  90  percent  with  an  average 
of  88.8  percent.  The  water  contents  measured  on  the  sediments 
adjacent  to  the  heater  probe  (sampled  by  the  over-core  and  tested  at 
SNL)  were  2  to  3  percent  below  the  water  contents  measured  in  the 
peripheral  cores.  These  lower  water  contents  may  be  due  to  drying 
during  the  delay  between  splitting  the  core  for  radiography  and 
subsampling  for  water  content.  The  water  content  at  the  mid-point 
of  the  heater  element  (0.25  m)  was  measured  as  6  to  7  percent  below 
the  water  content  values  at  0.20  and  0.30  m  (Figure  21).  This  low 
value  may  reflect  a  thermally-induced  consolidation  effect 
(Shephard,  1984). 

Radiographic  examination  of  the  over— core  was  performed  in 
order  to  look  for  sediment  cracking  resulting  from  heater  probe 
insertion.  The  results  of  this  examination  were  inconclusive. 

Cracks  developed  when  the  0.13  m  (5  in.)  diameter  over-core  was 
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split  making  it  very  difficult  to  differentiate  between  probe 
implant  cracks  and  core  splitting  cracks  (Shephard,  1984). 

The  undrained  shear  strength  measured  by  the  miniature  vane  was 
noted  to  be  reasonably  consistent  in  the  peripheral  core  tested  at 
SNL,  averaging  about  4.1  kPa  (Figure  22).  The  core  tested  at  NORDA, 
after  plane  and  auto  transport,  reflected  a  decrease  in  vane 
strength  with  depth  from  4.3  kPa  near  the  top  to  2.8  kPa  near  the 
bottom,  probably  reflecting  sediment  disturbance  due  to  transport. 
Sediment  sensitivities  as  measured  by  the  vane  are  4  and  greater. 

4.5  Strength  Prediction  from  Insertion  Pore  Pressures 


If  the  sediment  is  assumed  to  be  properly  modeled  as  an  elastic 
perfectly  plastic  material,  then  the  maximum  excess  pore  pressure  at 
the  probe  surface  is  given  by  Randolph  et  al.  (1979)  as: 


u 


imax 


s  In  (G/s  ) 
u  u 


(  14  ) 


Esrig  et  al.  (1977),  noting  that  the  ratio  G/su  can  often  be 
predicted  based  on  soil  type,  suggested  a  value  of  6  for  the  factor 
ln(G/ su)  in  Equation  (14)  ,  for  lean  inorganic  soils  of  moderate  to 
high  sensitivity,  resulting  in: 


In  ISIMU  II,  piezometer  probe  P-3  measured  on  insertion-induced 
excess  pore  pressure  of  11.6  kPa ,  and  vane  shear  strengths  measured 
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at  SNLA  on  a  peripheral  core  were  4.03  kPa  near  the  piezometer 
sensor  elevation.  The  resulting  pressure-to-strength  ratio  is: 

Uimax/Su  =  n‘6/4-03  “  2-9  *  3 . . 

This  low  value  of  the  uimax/su  ratio  is  measured  on  a  sample  of 
reconstituted  sediment,  and,  therefore,  it  is  not  considered  highly 
signif icant . 

5.0  SUMMARY  AND  CONCLUSIONS' 

The  surficial  sediment  crack  pattern  observed  in  the  ISIMU-II 
experiment  does  not  suggest  significant  penetration  of  the  cracks  to 
a  depth  where  the  excess  pore  pressure  and  heat  dissipation  might  be 
measureably  altered  by  the  cracking.  The  measured  similar  insertion 
and  dissipation  excess  pore  pressure  histories  for  the  two 
piezometers,  one  inserted  before  the  heater  probe  insertion,  the 
other  inserted  after,  lends  support  to  this  position. 

The  ISIMU-II  simulation  involved  the  collection  and  analysis  of 
data  pertaining  to  induced  excess  pore  water  pressure  build-up  and 
subsequent  dissipation  as  a  result  of  probe  insertion  and  thermal 
effects.  The  rate  of  dissipation  predicted  is  dependent  upon 
material  properties  as  well  as  the  dissipation  model  used.  The 
horizontal  coefficients  of  consolidation,  c^,  obtained  usincr 
Soderberg's  (1962)  method,  were  found  to  be  consistent  with  the 
ISIMU-I  results,  as  well  as  with  oedometer  test  results  (Silva  et 
al.,  1982).  The  c^  values  were  used  in  a  first  order  exponentially 
decreasing  dissipation  function  and  found  to  model  short  term 
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behavior.  However,  the  same  function  when  considered  in  terms  of 
long-term  performance  predicted  a  significantly  higher  dissipation 
rate.  The  predicted  behavior  was  "slowed  down"  using  Scott's  (1963) 
"smear"  factor  in  such  a  fashion  as  to  agree  with  a  least-squares 
approximation  to  the  data. 

The  thermally  induced  pore  pressures  were  modeled  using 
Mitchell's  (1976)  equation.  This  expression  involves  the  change  in 
pore  water  pressure  as  a  function  of  temperature  change  and  thermal 
properties  of  the  sediment.  The  actual  rate  of  temperature  change 
at  the  probe  was  not  known;  however,  these  data  were  available  for 
the  heating  element  and  container  walls.  Therefore,  the  general 
shape  of  the  curve  was  known  with  the  parameters  defining  the  curve 
assumed  a  priori.  The  thermally  induced  pore  pressure  model  was 
combined  with  the  dissipation  model  into  a  single  expression,  with 
the  former  serving  as  a  forcing  function  for  the  gain  in  pore 
pressure  and  the  latter  as  a  loss  term.  The  proposed  equation  was 
found  to  be  a  good  model  in  predicting  the  observed  thermally 
induced  pore  pressure  response. 

The  temperature  increase  calculated  from  the  temperature-excess 
pore  pressure  model  was  found  to  be  signif icantly  different  than 
that  which  would  have  been  predicted  from  ISIMU-I  results.  Possible 
causes  for  this  discrepancy  are  as  follows: 

a)  Improper  prediction  of  the  overall  behavior  of  the  sediment 
by  the  first  order  dissipation  model, 

b)  Excessive  "slowing  down"  of  the  dissipation  model  resulting 
■  in  a  lower  dissipation  rate. 
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c)  An  over  estimate  of  the  thermally  induced  pore  pressures  due 
to  errors  in  the  selection  of  thermal  coefficients, 

d)  Differing  heating  rates  and  temperatures  between  ISIMU-I  and 
ISIMU-II  simulations, 

e )  Differences  in  the  thermal  properties  and  response 
characteristics  of  the  pore  pressure  probes,  and 

f)  A  change  in  material  properties  as  a  function  of 

temperature,  namely,  changes  in  c,  . 

h 
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ADDENDUM  I.  Unadjusted  Data,  ISIMU-II  Experiment 
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27 

10 

23 

0 

23. 

1 

23. 

0 

23 

0 

23 

1 

23. 

5 

23. 

1 

-0 

1061 

0 

1849 

09 

27 

20 

23. 

1 

23. 

2 

23 

1 

22. 

0 

23. 

1 

23. 

5 

23 

1 

-0. 

1060 

0 

1810 

09 

27 

30 

23. 

1 

23 

1 

d  6  . 

1 

23. 

0 

23 

1 

23. 

5 

23 

1 

-0. 

1  059 

0. 

1771 

09 

27 

40 

23 

1 

23. 

1 

23. 

0 

23. 

1 

23 

2 

23. 

5 

23 

1 

-0 

1058 

0. 

1735 

09 

27 

50 

c.  6  . 

1 

23. 

CL 

23 

1 

23. 

0 

no 

d  w.«  . 

2 

23. 

5 

23 

1 

~0. 

1058 

0. 

1 6Q9 

09 

28 

00 

23. 

1 

23. 

2 

23. 

1 

22. 

0 

23. 

n 

d 

23. 

5 

23. 

i 

-0. 

1058 

0. 

1 66 1 

09 

28 

10 

23. 

0 

23. 

1 

23. 

1 

23. 

0 

23 

1 

23 

5 

23. 

Cl 

-0. 

1057 

0 

1627 

09 

28 

20 

23. 

0 

23. 

1 

n 

d  w>  . 

1 

nn 

d  w 

0 

23. 

1 

23. 

R 

23. 

1 

-0. 

1057 

0. 

1594 

09 

28 

30 

23. 

0 

23. 

1 

23 

1 

23 

1 

23 

1 

23. 

5 

'  23. 

1 

-0. 

1057 

0 

i  5 64 

09 

28 

40 

23. 

1 

23. 

2 

23. 

1 

23 

0 

23. 

1 

23. 

5 

23. 

1 

-0. 

1057 

0. 

1  534 

09 

28 

50 

23. 

0 

23 

1 

23. 

1 

22. 

0 

23. 

n 

c. 

23. 

5 

23. 

1 

-0. 

1058 

0. 

1  506 

09 

29 

00 

23. 

0 

23. 

1 

23. 

1 

23 

0 

23. 

1 

23. 

5 

23 

1 

-0. 

i  05E 

0. 

1478 

09 

29 

10 

23. 

0 

23 

1 

23. 

1 

23 

1 

23. 

1 

23 

5 

cl  6  . 

1 

-0. 

1059 

0. 

14  50 

09 

29 

20 

23. 

0 

23. 

2 

cl  6 

1 

no 

d  o . 

0 

23 

1 

23. 

5 

23 

1 

-0 

1059 

0 

1424 

09 

29 

30 

23 

1 

23. 

1 

23. 

1 

n  o 

c.  6  . 

0 

23. 

2 

23. 

5 

23 

1 

-0. 

1059 

0. 

1399 

09 

29 

40 

23, 

0 

23 

1 

23. 

0 

23 

0 

23. 

2 

23. 

5 

23. 

1 

-0. 

1059 

0 

1374 

09 

29 

50 

23. 

0 

23 

1 

23. 

1 

23 

0 

23 

cl 

23 

5 

23. 

1 

-0. 

1060 

0 

1349 

09 

30 

00 

23 

0 

23. 

2 

23. 

1 

23. 

0 

23 

1 

23. 

5 

23. 

i 

“0. 

1 060 

0 

1  327 

09 

30 

10 

23. 

1 

23. 

cl 

23. 

1 

23. 

1 

23 

2 

23 

5 

no 

d  w* 

d 

-0. 

1061 

0. 

1304 

09 

30 

20 

23 

1 

23. 

1 

23. 

1 

23. 

0 

23. 

i 

23. 

r 

23. 

l 

-0. 

1061 

0 

i  282 

09 

30 

30 

23. 

0 

23. 

1 

23. 

1 

23 

0 

23. 

cl 

23. 

6 

no 

d  w<  . 

1 

-0 

1061 

0. 

1261 

09 

30 

40 

23 

0 

23. 

2 

23 

0 

23. 

0 

23. 

2 

no 

d  d  . 

5 

23. 

l 

-0 

1  06d 

0 

1239 

09 

30 

50 

23. 

1 

23. 

1 

23. 

0 

23. 

0 

23. 

2 

23. 

5 

23 

i 

““0. 

1063 

o 

1  d  i v 

09 

31 

00 

23. 

0 

23. 

1 

23. 

0 

23 

0 

23. 

1 

23. 

5 

23. 

i 

-0. 

1062 

0 

119- 

09 

31 

10 

23 

0 

23 

1 

23. 

1 

23. 

0 

23 

1 

27 

5 

23. 

l 

-0 

1063 

A 

1  180 

09 

31 

20 

23. 

1 

23. 

1 

23. 

0 

23. 

0 

23. 

cl 

23 

4 

23. 

l 

“0. 

1064 

0 

1  IcO 

09 

31 

30 

23. 

1 

23. 

2 

23. 

1 

23. 

1 

23. 

n 

d 

n  o 

d  w  . 

e. 

•u 

23. 

l 

-0. 

1064 

0 

1141 

09 

31 

40 

23. 

1 

23 

n 

c. 

23. 

1 

23. 

1 

23 

O 

d 

23. 

5 

23 

l 

-0. 

1064 

0 

1126 

09 

31 

50 

23. 

1 

23. 

n 

23. 

1 

23 

0 

23. 

1 

23. 

4 

23. 

e_ 

“0. 

1064 

0. 

1  105 

09 

32 

00 

23 

0 

23. 

1 

23 

0 

23. 

0 

23 

1 

no 

d%j 

4 

23. 

1 

-0 

1064 

0 

108"7 

09 

32 

10 

23. 

1 

23. 

2 

23. 

1 

23. 

0 

23. 

1 

23. 

4 

23. 

1 

-0. 

1065 

0. 

1071 

09 

32 

20 

23. 

1 

23. 

2 

23. 

1 

23 

0 

23. 

1 

23, 

4 

no 

d  —> 

4 

X 

-0. 

1065 

0 

1054 

09 

32 

30 

23 

1 

23. 

2 

23 

1 

23. 

0 

23. 

O 

c. 

23. 

4 

23. 

1 

“0. 

1066 

0 

1038 

09 

32 

40 

23. 

1 

23 

2 

23. 

0 

23. 

1 

23. 

2 

23. 

3 

23. 

1 

-0. 

1066 

0. 

,  1021 

09 

32 

50 

23 

1 

23. 

2 

23. 

0 

23. 

0 

23. 

o 

CL 

23. 

3 

23. 

1 

-0. 

1066 

0 

1005 

09 

33 

00 

23. 

1 

23 

2 

23. 

1 

23. 

0 

23 

1 

23 

3 

23. 

1 

-0. 

1066 

0. 

0990 

09 

33 

10 

23. 

1 

23. 

2 

23. 

1 

23 

0 

23. 

1 

23. 

n 

c. 

23 

1 

-0 

1066 

0 

097  5 

09 

33 

20 

23. 

1 

23. 

2 

23. 

1 

23. 

0 

23 

1 

23. 

2 

23. 

1 

“0 

1067 

0 

.  0959 

09 

33 

30 

23 

1 

23 

2 

23. 

1 

23 

0 

23 

1 

23. 

2 

23. 

1 

-0. 

1067 

0 

0944 

09 

33 

40 

23 

1 

23. 

2 

23. 

1 

23. 

0 

23. 

2 

23. 

1 

23 

1 

-0. 

1067 

0 

.  0930 

09 

33 

50 

23 

0 

23 

2 

d  vd  . 

1 

23. 

0 

23 

1 

23 

i 

dJ . 

1 

-0. 

1067 

0 

.  09x7 

09 

34 

00 

23. 

1 

23. 

n 

c. 

23. 

0 

23 

0 

23 

2 

23. 

l 

23. 

1 

-0 

1068 

0, 

.  0903 

09 

34 

10 

23. 

1 

23. 

2 

23. 

1 

23. 

0 

23. 

2 

23 

0 

23. 

1 

-0. 

1066 

0 

08E9 

09 

34 

20 

23 

1 

23. 

1 

23 

0 

23. 

0 

23. 

1 

23. 

.  0 

23. 

1 

-0 

1068 

0 

.  087 o 

09 

34 

30 

23. 

0 

23. 

1 

23. 

1 

23 

0 

23 

1 

22. 

9 

23. 

1 

“0. 

1069 

0 

0861 

09 

34 

40 

23. 

0 

23. 

1 

23. 

1 

23. 

0 

23 

i 

22. 

9 

23. 

1 

“0. 

1069 

0 

.  0848 

09 

34 

50 

23. 

1 

23. 

2 

23. 

1 

23 

0 

23, 

2 

23, 

,  0 

23. 

n 

d 

-0. 

1069 

0. 

.  0835 

09 

35 

00 

23. 

1 

23 

1 

23. 

1 

23. 

0 

23. 

1 

22. 

.  8 

23. 

1 

-“0. 

1069 

0. 

.  0823 

09 

35 

10 

23. 

0 

23. 

1 

23 

1 

23 

0 

23. 

1 

22. 

.  6 

23. 

1 

“0. 

1069 

0 

.  0810 

Thermistors 


Hours 

Min 

Sec 

0 

1 

2 

3 

4 

5 

6 

P-3 

P-2 

09 

35 

20 

23,  1 

23.  2 

23.  1 

23.  0 

23.  1 

22.  S 

23.  1 

-0  1069 

0. Q7Qo 

09 

35 

30 

23.  1 

23  2 

23  1 

23.  0 

23.  2 

22.  8 

23  1 

-0. 1070 

0  0785 

09 

35 

40 

23.  1 

23.  1 

23.  0 

23  0 

23.  2 

22.  7 

23.  1 

“0.  1070 

0.  0774 

09 

35 

50 

23.  1 

23  2 

23  1 

23  1 

23  2 

22.  8 

23.  1 

-0.  1070 

0  0762 

09 

36 

00 

23.  1 

23.  1 

23.  0 

23.  0 

23  1 

22.  7 

23  1 

-0.  1070 

0.  0751 

09 

36 

10 

23.  1 

23.  2 

23.  1 

23  0 

23.  2 

22.  7 

23.  2 

-0  1071 

0  074  0 

09 

40 

00 

23.  1 

23  2 

23.  0 

23,  0 

23.  1 

22.  2 

23.  1 

-0.  1074 

0  0513 

09 

45 

00 

23.  0 

23.  1 

23  0 

23.  0 

23.  1 

22  3 

23  1 

-0.  1078 

0.  0290 

09 

50 

00 

23.  0 

23.  1 

23.  0 

22.  9 

23.  1 

23  0 

23.  1 

-0,  1086 

0  0113 

09 

55 

00 

23.  1 

23.  1 

23  0 

23  0 

23.  1 

23.  2 

23  1 

-0.  1096 

-0  0031 

10 

00 

00 

23.  0 

23  1 

23  0 

23.  0 

23,  1 

23.  3 

23.  1 

-0.  1109 

-0  0151 

10 

05 

00 

23  1 

23.  1 

23  0 

23  0 

23.  2 

23.  4 

23.  1 

-0  1121 

-0  0250 

10 

10 

00 

23.  1 

23  2 

23  0 

23  0 

23  2 

23.  4 

23.  1 

-0,  1135 

-0  0338 

10 

15 

00 

23.  0 

23.  1 

23  1 

23  0 

23.  1 

23  5 

23.  1 

-0.  1149 

-0  041 1 

10 

20 

00 

23.  0 

23.  1 

23.  0 

23  0 

23  1 

23  1 

23.  1 

-0.  1  161 

-0  0474 

10 

25 

00 

23.  0 

23  1 

23  0 

23  0 

23  1 

22.  4 

23.  1 

-0.  1173 

-0.  0532 

10 

30 

00 

23.  0 

23.  1 

23  0 

23.  0 

23  1 

22.  1 

23.  1 

-0.  1  185 

-0.  0584 

10 

35 

00 

23  0 

23,  1 

23.  0 

23.  0 

23  2 

22.  5 

23.  1 

-0.  1192 

-0.  0630 

10 

40 

00 

23.  0 

23  1 

23.  0 

23  0 

23.  1 

22  9 

Cm3  1 

-0  1203 

-0.  067 1 

10 

45 

00 

23  0 

23  2 

23  1 

23.  0 

23.  2 

23.  2 

23  1 

-0  1213 

-0.  0708 

10 

50 

00 

23.  1 

•  23.  2 

23  1 

23  1 

23.  2 

23.  3 

23  1 

-0  1223 

-0  0741 

10 

55 

00 

23.  0 

23  1 

23.  1 

23  0 

23.  1 

23.  4 

23.  1 

-0.  1234 

-0  0771 

1 1 

00 

00 

23  0 

23.  1 

23  1 

23  0 

23  1 

23.  4 

23  1 

-0  1244 

-0  0797 

1 1 

05 

00 

23.  0 

23  1 

22  9 

22  9 

23  1 

23  5 

23  3 

-0  1253 

-0  0818 

1 1 

10 

00 

23.  1 

23  1 

23.  1 

23.  0 

23  1 

23.  7 

22  1 

“0.  1263 

-0.  0845 

1 1 

15 

00 

23  1 

23.  1 

23.  1 

23.  0 

23.  1 

23.  1 

23  1 

“0  1272 

-0. 0866 

11 

20 

00 

23.  0 

23.  1 

23.  0 

23  0 

23.  1 

22.  7 

23.  1 

-0.  1280 

“0  0S8c> 

1 1 

20 

30 

23.  1 

23  2 

23.  1 

23  0 

23.  1 

22.  7 

23.  2 

-0.  1281 

—0,  0357 

1 1 

21 

30 

23.  0 

23  1 

23  0 

23  0 

23.  1 

22.  5 

23.  1 

-0. 12S3 

-0.  0890 

11 

22 

30 

23.  1 

23  2 

23  0 

23.  0 

23.  1 

22.  6 

23.  1 

-0.  1284 

— 0.  0394 

DATA 

SET 

4: 

HEATER  PROBE  TURNED  ON 

11 

23 

00 

23.  3 

23  3 

22  7 

22.  8 

23  4 

22  7 

23.  0 

-0.  1285 

-0  08-2 

1  1 

23 

30 

23.  4 

24  4 

24.  S 

25  1 

25  2 

22  4 

22  8 

-0  1284 

-0  089- 

1  1 

24 

00 

23.  1 

25.  6 

27.  2 

27.  2 

27.  0 

22  2 

23.  0 

-0.  1288 

“0.  0398 

1  1 

24 

30 

23.  4 

27.  2 

29.  0 

28.  5 

29.  0 

22.  8 

23.  0 

-0.  1296 

-0.  09C 3 

11 

25 

00 

23.  2 

27.  5 

30.  1 

29.  7 

30.  1 

22.  1 

22.  8 

-0.  1291 

-0.  0906 

11 

25 

30 

23.  3 

28.  9 

31  5 

30.  7 

31  5 

22.  5 

22.  9 

-0.  1299 

“0. 0909 

1  1 

26 

00 

23  3 

29  7 

32  6 

31  3 

32.  6 

22.  7 

22.  9 

-0.  1296 

-0.  09 1 0 

11 

27 

00 

23.  5 

30.  8 

34.  5 

33.  4 

34  8 

2.d  3 

22.  8 

“0.  1290 

-0  0914 

11 

28 

00 

23.  5 

32  0 

36  2 

34.  9 

36.  4 

22.  3 

22.  8 

-0.  1282 

“0.  0913 

1  1 

29 

00 

23  1 

32  8 

37.  8 

36.  1 

37.  5 

22.  0 

23.  0 

-0  1273 

-0.  09 OS 

11 

30 

00 

23  2 

33  7 

38  9 

37.  2 

38.  7 

22  0 

22.  9 

-0.  1262 

-0  0907 

11 

31 

00 

22  8 

34  6 

40.  6 

38.  1 

39.  5 

n 

Cm  Cm  .  Cm 

23.  3 

-0.  1253 

— 0.  0905 

1  1 

32 

00 

22  8 

35  4 

41  6 

39.  0 

40.  6 

22  2 

23  1 

-0  1241 

-0  0900 

1  1 

33 

00 

23.  2 

36.  3 

42.  1 

39  5 

42  0 

22.  4 

22.  9 

-0  1234 

-0.  08°  1 

1  1 

34 

00 

23.  0 

36.  5 

43.  1 

40  7 

42.  6 

22.  0 

23.  1 

-0. 1222 

-0, 088- 

1  1 

35 

00 

22.  7 

37.  4 

44  0 

41.0 

43.  0 

22  4 

23.  3 

-0.  1209 

-0  0883 

1  1 

36 

00 

22.  6 

38.  0 

44.  9 

41.  5 

43  6 

22.  6 

23.  3 

-0.  1203 

-0  0877 

11 

37 

00 

23.  3 

38  9 

45.  1 

42.  2 

45.  1 

22.  9 

22.  9 

-0.  1  194 

-0.  0870 

11 

39 

00 

22.  7 

39  5 

46.  8 

43.  4 

45.  8 

23.  0 

23.  3 

-0. 1175 

-0. 0863 

11 

41 

00 

22.  8 

40.  2 

47.  9 

44.  6 

47.  0 

22.  9 

23  c: 

-0.  1163 

“0.  0854 

1  1 

43 

00 

22.  8 

41.  3 

48  9 

45.  2 

47.  9 

23.  3 

23.  2 

-0  1147 

— 0  0845 

54 

57 

07 

19 

31 

39 

49 

59 

14 

24 

34 

44 

54 

04 

14 

24 

34 

44 

54 

04 

14 

24 

34 

10 

20 

30 

40 

50 


Thermistors 


Sec 

0 

i 

2 

3 

4 

5 

6 

P-3 

00 

22  8 

42.  5 

50.  1 

45  9 

49.  0 

23  9 

23.  5 

“0.  1141 

00 

23.  5 

43.  1 

50.  3 

46  7 

50.  4 

23.  8 

23.  0 

“0.  1132 

00 

23.  4 

44.  1 

51.  5 

47.  7 

51.  5 

23.  9 

23  1 

-0  1118 

00 

22.  9 

45  2 

53.  3 

48.  7 

52.  2 

24.  1 

23 .  6 

-0. 1101 

00 

23.  1 

45.  2 

53  6 

49.  7 

53.  2 

23.  4 

22  9 

— 0  1089 

00 

23.  1 

48.  0 

56.  2 

51  3 

55.  6 

24  3 

23  3 

-0.  10fc9 

00 

23.  1 

49.  2 

58  1 

53  8 

57.  8 

22.  3 

23.  1 

-0  1051 

00 

22  9 

50.  7 

59.  7 

55  1 

59.  4 

22.  0 

23.  0 

-0  1036 

00 

23  0 

52.  3 

60.  8 

55  5 

60.  5 

22.  7 

c-3  1 

-0. 1033 

00 

22  8 

52  3 

62  0 

57.  1 

61.  1 

22.  5 

23.  3 

-0. 1028 

00 

23  0 

53  9 

63.  0 

57.  5 

62.  3 

23.  4 

23.  3 

“0  1033 

00 

23  0 

54  3 

63.  6 

58  7 

63  4 

22  8 

23.  0 

— 0.  1039 

00 

23  0 

55.  5 

64.  7 

59  1 

64  1 

23.  5 

23.  4 

-0.  1052 

00 

23.  0 

56.  1 

65.  4 

59  7 

64  7 

23.  7 

23.  5 

“0.  1066 

00 

23.  0 

56.  0 

65.  4 
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ADDENDUM  II.  Reduced  Data,  ISIMU-II  Experiment 
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ADDENDUM  III.  Excess  Pore  Pressure  Prediction  Program  MENUDO 


Hi  SB!  ARE  CURVE  h IT  PARAMETERS  FOR  TEMP . -f ( TIME)  RELATIONSHIP 


’  U  !  -’RUGRAM  "MENUDO” 

20  !  THIS  PROGRAM  PLOTS  CALCULATES  EXPECTANT  PORE  PRESSURES  BASED 
: IME  lFFlCTS 
20  OPTION  BASE  ! 

AO  !  PRINTER  IS  701.80 

^0  DIM  M  <  50 )  .T<50)  .Li!  (50)  .02(50)  .03(50)  .04(50)  .X(50)  ,  Yl  (50)  ,  Y2<50> 
60  DIM  D ( 50 ) . T4 ( 50 ) . U( 50 ) . P ( 50 ) 

70  M=8 
80  N=47 
90  ! 

100  A 1 =  2 1 0 
110  B . = . 25 

i 20  !  70-INITIAL  TEMPERATURE 
130  70=22.8 

i 40  !  N ! = INITIAL  POROSITY 
1 50  N1=.74 

160  !  A2=7HERMAL  COEFFICIENT  OF  EXPANSION  FOR  SOIL  SOLIDS 
170  A2=. 000001 

:80  !  A3= THERMAL  COEFFICIENT  OF  EXPANSION  FOR  SOIL  STRUCTURE 
190  A3--. 00005  PER  °C_ 

200  !  M ; =  Mvs  COEr FICIENT  OF  VERTICAL  COMPRESSIBILITY  FOR  REBOUND 
210  Ml =-.00025 

220  !  C !  =Cv/  COEFFICIENT  OF  CONSOLIDATION/HYDRAULIC  DIFFUSIVITY 


230 

!  FOR  PROBE  P-3  IN 

CM' 2/M IN 

240 

Cl =.01 9277 

250 

!  FOR  PROBE  P-2  IN 

CM  '  2/MIN 

260 

!  Cl =.014035 

270 

!  R1 = PROBE  RADIUS 

280 

Ri  =  .4 

290 

!  M2=SMEAR  FACTOR/I 

NTERFACE 

300 

!  FOR  PROBE  P-3 

310 

M2 = 4.6 

320 

!  FOR  PROBE  P-2 

330 

1  fvj  -  C*  1> 

340  !  PORE  PRESSURE  PARAMETERS  DUE  TO  PREVIOUS  HISTORY 

350  !  FOR  PROBE  P-3 

360  A4=  .  !  898  •«  B4  =  -.006l93  i>  T4=1Q0 

370  !  FOR  PROBE  P-2 

380  !  A4=.5856  w  B4  =  -.01726  14=100 

390  C4=A4*EXP(B4*-T4) 

400  !  READS  IN  PREVIOUSLY  STORED  DATA 
410  FOR  Q=M  TO  N 

420  ASSIGN*  !  TO  "PR0BEDATA4 : D70 1 " 

430  READ-  1.0  :  M(Q) . Tl . T2 . T3 . T4<Q> . T5 . T6 . T7 , U1 (0 ) . U2< 0 ) 
440  ASSIGN*  1  TO  * 

450  NEXT  Q 
460  GOTO  530 

470  DISP  "DO  YOU  WANT  TO  CHANGE  AdB  PARAMETERS  Y/N"  *  BEEP 
480  INPUT  BS 


490  IF  BS="N"  THEN  1200 

500  DISP  "INPUT  AdB"  >1  BEEP 

510  INPUT  A1 .81 

520  PRINT  A 1 ; B 1 

530  0=0 

540  U=0 

550  FOR  I=M  TO  N 
560  0=0+1 


570  X(0)=M( I)-M(M) 

580  !  PROBE  P-3  VOLTAGE-Tu-PRESSURE 
590  Y 1 <  0 )  =  . 0 1 1 1  +  ( U1 < I )  +  . 1 505  >*3 . 998 
800  !  PROBE  P-2  VOLTAGE-TO-PRESSURE 


CONVERSION 

CONVERSION 


ON 


610  !  Y1  <0>-.0391+<U2<I>+.107 1  >*3.995 
620  !  TEMPERATURE  AS  A  FUNCTION  OF  TIME 
630  T(O)*X(O)/(A1+Bl«X(O))+T0 

640  1  AO=THERMAL  COEFFICIENT  OF  EXPANSION  FOR  SEA  HATER 
650  A0= .000064609+. 00001 0474*T<0)- . 000000056*1 ( 0> ' 2 
660  IF  0=1  THEN  T2  =  0  ELSE  T2=T(0>-T<0- 1 ) 

670  !  U3=P0RE  PRESSURE  INCREASE  DUE  TO  THERMAL  EFFECTS 
680  U3(0)  =  T2/M  WN1*<A2-A0)+A3> 

680  !  U4=P0RE  PRESSURE  INCREASE  CORRECTION  DUE  TO  PREVIOUS  HISTORY 
700  U4(0)=C4-A4*EXP(B4*(T4+X(0) ) > 

710  U=U3(0>+U 
720  U(0>=U+U4(0> 

730  !  PI =PORE  PRESSURE  DISSIPATION  AS  A  FUNCTION  OF  TIME 
740  B0=-(C1/R1'2/M2> 

750  P=0 

760  FOR  L=1  TO  0 

770  IF  L= 1  THEN  PI -0  ELSE  P 1 =U3< L- 1 )*EXP(B0* < X< 0 ) -X< L- 1 ) ) ) 

780  P-P+P1 
790  NEXT  L 

800  IF  0=1  THEN  P(O)-YKI)  ELSE  P( 0 > =U3< 0> +P+ Y 1 < 1 > 

810  IMAGE  DD , 2X . DDD . DO , 2X . D . DDD  .2X.  D.DDD.2X.D.DDD 

820  !  PRINT  USING  810  :  0 , X ( 0  > . Y 1 <  0  > .  U  ( 0  >  +  Y 1 < 1 > . P <  0  >  +  Y 1 < 1 > 

83C  IMAGE  D.DDDD.2X.D.DBBD.2X.B .DDD 
840  !  PRINT  USING  830  ;  U3<0> .U4(0) ,P<0> 

850  !  PRINT 
860  NEXT  I 
870  DEG 

880  PLOTTER  IS  705 
890  LIMIT  10.230.10.180 
900  LOCATE  10.125,10.95  . 

910  SCALE  0.300.0.1 .2 
920  FXD  0.1 

930  LAXES  -25. . 1 ,0.0. 2, 2 
940  CSIZE  4 
950  SETGU 

960  MOVE  68.5  lORG  5 
970  LABEL  ’’TIME  (MIN.)" 

980  MOVE  3.53  lDIR  90 
990  LABEL  "PRESSURE  <P$I>" 

1000  SETUU  »  LORG  2  ¥  LDIR  0 

1010  MOVE  150,1.2  €  LABEL  "PORE  PRESSURE  PROBE  P-3" 

1020  !  MOVE  150.1.2  *¥  LABEL  "PORE  PRESSURE  PROBE  P-2" 

1030  MOVE  150.1.15  «  LABEL  "*  -  DATA" 

1040  MOVE  150.1.1  $  LABEL  "+  -  THERMALLY  INDUCED  PORE  PRESSURE" 
1050  MOVE  150,1.05  *  LABEL  "0  -  THEORY  PREDICTION" 

1060  LORG  5  C«  CSIZE  2 

1070  FOR  J=1  TO  0 

1080  MOVE  X< J) . Y1 < J)+U4( J> 

1090  LABEL  "*" 

1100  NEXT  J 

1110  FOR  K=1  TO  0 

1120  MOVE  X<K)  ,U<K)-U4(K)+YH1 ) 

1130  LABEL  "+" 

1140  NEXT  K 

1150  FOR  L=1  TO  0 

1160  MOVE  X(L  > . P ( L ) 

'1170  LABEL  "0" 

1180  NEXT  L 
1190  GOTO  470 

12CC  DISP  "ADIDS  SENDR  RIGONES" 

1210  END 


Figurel.  ARRANGEMENT  FOR  SEDIMENT  CONSOLIDATION 

AT  SNLA 


Figure  2.  PLAN  VIEW  OF  TEST  ARRANGEMENT  (SCALE  1  :4) 
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Figure  3.  PIEZOMETER  PORE  PRESSURE  RESPONSE  FOR  THE  COMPLETE  TEST  PERIOD 

(  lpsi  =  6.895kPa  ) 
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Figure  4.  PIEZOMETER  PORE  PRESSURE  RESPONSE  DURING  EVENT  1:  PROBE  P-3  INSERTION 
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Figure  5.  PIEZOMETER  PORE  PRESSURE  RESPONSE  AT  PROBE  P-3  DURING  EVENT  2:  HEATER  INSERTION 

(  lpsi  =  6.895kPa  ) 
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Figure  6.  PIEZOMETER  PORE  PRESSURE  RESPONSE  DURING  EVENT  3:  PROBE  P 
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Figure  7.  PIEZOMETER  PORE  PRESSURE  RESPONSE  DURING  EVENT  4:  HEATER  ACTIVATION 

(  lpsi  =  6.895kPa  ) 
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Figure  8.  GRAPHICAL  REPRESENTATION  OF 
PORE  PRESSURE  GAIN-LOSS  MODELS 
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Figure  9. 
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THE  LOG  (TIME)  DISSIPATION  CURVE  FOR  PROBE  P-3 
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Figure  10.  THE  LOG  (TIME)  DISSIPATION  CURVE  FOR  PROBE  P-2 
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Figure  11.  THE  LOG  (TIME)  DISSIPATION  CURVE  AT  PROBE  P-3  DUE  TO  THE  INSERTION  OF  HEATER  PROBE 


Figure  12.  THE  DISSIPATION  CURVE  AT  PROBE  P-3 
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Figure  13.  THE  DISSIPATION  CURVE  AT  PROBE  P-2 
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Figure  14.  THE  PREDICTED  DISSIPATION  CURVES  AT  PROBE  P-3  AS  A  FUNCTION  OF  THE  SMEAR 
FACTOR ,  m 
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PORE  PRESSURE  PROBE  P-2 
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Figure  15.  THE  PREDICTED  DISSIPATION  CURVES  AT  PROBE  P-2  AS  A  FUNCTION  OF  THE  SMEAR 
FACTOR,  m 
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Figure  16.  HISTORY  CORRECTION  FOR  THE  THERMALLY  INDUCED  PORE  PRESSURES  AT  PROBE  P-3 
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Figure  17.  HISTORY  CORRECTION  FOR  THE  THERMALLY  INDUCED  PORE  PRESSURES  AT  PROBE  P-2 

(  lpsi  =  6.895kPa  ) 
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Figure  18.  TEMPERATURE  AT  THE  HEATING  ELEMENT  AS  A  FUNCTION  OF  TIME 
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Figure  19.  MODEL  PREDICTIONS  AT  PROBE  P-3  (  lpsi  =  6.895kPa  ) 
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Figure  20.  MODEL  PREDICTIONS  AT  PROBE  P-2  (  lpsi  =  6.895kPa  ) 
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Figure  21.  WATER  CONTENTS  OF  TEST  TANK  SEDIMENTS 
AT  CONCLUSION  OF  ISIMU-II  EXPERIMENT 
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Figure  22.  UNDRAINED  SHEAR  STRENGTHS  OF  TEST  TANK  SEDIMENTS 

AT  CONCLUSION  OF  ISIMU-II  EXPERIMENT 
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